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Abstract 


The effect of inhomogeneous compression of GDL on the mass and charge transfer in PEMFC is studied. The model utilizes experimentally 
evaluated GDL parameters as a function of thickness. The modeling results are compared with a conventional model that excludes the effects. As 
a result, it is shown that the inhomogeneous compression has a significant effect on the current density distribution because of the varying contact 
resistance between GDL and electrode. This also implies that there are possible hot spots occurring inside the electrode, and thus inhomogeneous 
compression can have significant effects on the lifetime and local performance of the cell. According to the achieved results, the inhomogeneous 


compression of GDL cannot be neglected. 
© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Fuel cells are electrochemical devices that convert the chem- 
ical energy of reactants directly into electricity and heat. Due 
to their advantageous properties, such as potential for high 
energy density and low environmental emissions, fuel cells 
are believed to gain significant market in the near future. The 
main applications for fuel cells are automotive, stationary, and 
portable power production. The large-scale market penetration 
of fuel cells still requires cost and performance improvements. 
In order to achieve these improvements, it is essential to have a 
deep insight into the processes occurring inside the cell and its 
components. 

One of the key components affecting the performance of a 
polymer electrolyte membrane fuel cell (PEMFC) is the gas dif- 
fusion layer (GDL). GDLs have to provide several functions for 
the fuel cell operation: a passage for reactant access and excess 
product water removal to and from the electrodes, electronic 
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conductivity, heat removal, and adequate mechanical support 
for the membrane electrode assembly (MEA). GDLs are typi- 
cally made of highly porous carbon-fiber based paper or cloth in 
order to fulfill these requirements. High porosity gives a char- 
acteristic soft and brittle structure for the GDLs, which causes 
a deformation in its shape when the fuel cell is assembled and 
components compressed together. 

The physical properties of GDL are changed under compres- 
sion, and thus also its mass, heat, and charge transfer properties 
are changed. Any change in physical properties of GDL in 
order to improve the charge transport may cause an adverse 
effect on the mass transport and vice versa. It has been exper- 
imentally shown that changes in the properties can have a 
significant effect on the fuel cell performance, see e.g. [1,2]. 
It is particularly worth noting that the deformation of GDL is 
not homogeneous. The parts of the GDL situated under the 
current collecting rib of the flow-field plate are significantly 
more compressed than the parts under the channel. This inho- 
mogeneous compression may cause significant changes in the 
local physical properties of GDL, and thus also in the local 
cell performance by changing the local reactant and current 
profiles. 

Even though it has been known for some time that the inho- 
mogeneous compression of GDL may have a significant effect 
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Nomenclature 


N bh CLIN 2S BAD mero Sy th af 


concentration (mol m73) 
channel and rib width (m) 
diffusion coefficient (m? s7!) 
unit vector 

reversible cell potential (V) 
Faraday constant (A s mol!) 
thickness (m) 

current density (A m~?) 
current production rate (A m~?) 
exchange current density (A m7?) 
permeability (m7) 

molar mass (kg mol” !) 

molar flux (mol m~? s7!) 
pressure (Pa) 

resistance (Q m2) 

gas constant (J mol`! K7!) 
liquid water saturation 
temperature (K) 

velocity (m s7!) 

molar fraction 

number of transferred electrons 


Greek letters 


Or reaction symmetry factor 
E porosity 

n overpotential (V) 

u viscosity (kg m7! s7!) 

p density (kg m~?) 

o conductivity (Q7! m7!) 
O electronic potential (V) 
dm protonic potential (V) 
Subscripts 

a anode 

ave average 

c cathode 

comp compressed 

cont,e contact between GDL and electrode 
cont,gr contact between GDL and graphite 
e electrode 

eff effective 

GDL gas diffusion layer 

H20 water 

min minimum 

N2 nitrogen 

O2 oxygen 

ref reference 

sat saturation 

x x-direction, in-plane 

y y-direction, through-plane 


on the cell performance, most of the PEMFC modeling stud- 
ies have neglected this effect. Zhou et al. [3] investigated the 
effect that the shape of the current collecting rib has on the 
changes in porosity of GDL and contact resistance between rib 
and GDL. Sun et al. [4] studied the effect that the inhomoge- 
neous compression of GDL, affecting the local porosity and 
conductivity, has on the fuel cell performance and local current 
density distribution. Even though they assumed a fairly small 
value for the compression (15%) and neglected the permeability 
and contact resistance effects, they concluded that the inhomoge- 
neous compression affected the local current density distribution 
notably. Sui and Djilali [5] varied only the values of through- 
plane conductivity and diffusivity of GDL, and observed 
that the changes in local values affect the current density 
profile. 

This second part of this contribution focuses on modeling 
the effect that inhomogeneous compression of GDL has on local 
species and current distributions. The model utilizes experimen- 
tally evaluated parameter values as a function of GDL thickness. 
These values are taken from the first part of this study, which 
focused on the ex-situ experimental evaluation of the GDL 
parameters [6]. The modeling results are compared with a con- 
ventional model that excludes the effects of inhomogeneous 
compression, and assumes the GDL parameters constant. The 
comparison gives insight into how the inhomogeneous compres- 
sion of GDL affects the local cell performance. 


2. Model 


Two different cases are modeled: one with homogeneous 
properties of GDL (referred to as “base case’) and one where 
the inhomogeneous compression of GDL is taken into account. 
The used geometry is a 2D cross-section of the cell, and the mod- 
eled geometries are illustrated in Fig. 1. The model consists of 
the anode and cathode GDLs and electrodes, and the membrane. 
The ribs and channels of the flow-field plates are accounted for 
as boundary conditions. Only half-widths of the rib and chan- 
nel structure and components below them are modeled, and the 
left and right geometry edges of Fig. 1 (boundaries IN-VII) 
are modeled with symmetry boundary conditions, i.e. it is 
assumed that the cell geometry continues symmetrically to both 
directions. 

The model takes into account the charge and multicomponent 
mass transfer in the cathode GDL and electrode, and charge 
transfer in the membrane and anode GDL and electrode. The 
main assumptions of the model are that water may exist in 
two phases but the transfer of liquid phase is similar to gas 
phase, i.e. equations for capillary movement are not included. 
In addition, the anode activation and mass transfer limitations 
are assumed to be negligible, and the cell is treated as isother- 
mal. Even though the effects of inhomogeneous compression are 
taken into account also at the anode, the intrusion of the GDL 
into the channel is not included in the modeled geometry. This is 
made for simplicity, because the inclusion of it has an insignif- 
icantly small effect on the current profile of the anode GDL 
only. The details of the model are described in the following 
subchapters. 
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Fig. 1. Modeling domains: (a) geometry with homogeneous compression of GDL (base case); (b) geometry with inhomogeneous compression of GDL taken into 
account. Roman numerals refer to the boundaries of the modeled geometries. The domains representing different cell components are not in scale in the figure. 


2.1. Equations 


2.1.1. Cathode GDL 

The governing equations at the cathode GDL are the conser- 
vation of mass, momentum, species, and charge listed in Eqs. 
(1)-(4), respectively: 


V- (pv) =0 (1) 
U- 
Vp=-—— 2 
P r” 2) 
V-Ñi=0 (3) 
IPGDL,c = O@GDL,c = 
V + | —OGDL,x Èy — OGDL,y <2, )=0 (4) 
Ox oy 


where different in-plane and through-plane conductivities of 
GDL are presented with subscripts x and y, respectively. 

The multicomponent mass transfer of different species (oxy- 
gen, water, and nitrogen) takes into account the convective and 
diffusive mass fluxes. The species flux equation is 


No, = || Xo = | VXo, 
> = CDeff (5) 
N50 X50 V Xm o 


Concentration and density of the gas mixture are calculated 
from the ideal gas law: 


P 
c= RT (6) 
and 
pea isl (7) 
RT 
where M is the molar mass of the gas mixture defined as 
M =J XM; (8) 
i 


Deg in Eq. (5) is the effective multicomponent diffusion coef- 
ficient tensor corrected by the Bruggeman correlation to take the 
effect of porosity and tortuosity into account: 
Dest = (e(1 — s))'°D (9) 

The effect of porosity reduction due to liquid water satu- 
ration, s, is also accounted for. The saturation s is defined as 
the fraction of pores occupied by liquid water, i.e. the molar 
fraction of water exceeding the corresponding saturation molar 
fraction: 


Psat 


Xs = (10) 


The saturation pressure of water can be calculated as [7]: 


log o(Psat (bar)) = 28.59051 — 8.2 log(T + 0.01) 


3142.31 


0.0024804(T + 0.01) — 2 
i U t00- F001 


(1) 
which gives the saturation pressure in bar. The components 
of multicomponent diffusion coefficient tensor D are calcu- 
lated from binary Maxwell—Stefan diffusion coefficients as 


[8]: 


Xo, Dy30,N) + (1 — X0,)D0>,H20 


Di = Doy,N> 5 
Do,,N, — Dos,H0 
D2 = Xo, Damo, N, —— 5 =e as 
Dy50,N. — Doo, 
Diu = Xmo Do, n ——— F a, 
X4,0Do,,N, + (1 — Xm o)Do:,m 0 
Dz = Dy,0,N> 5 ; 


S = Xo, Dy,0,N, + Xm o Do, Na + XN, Do, m0 (12) 
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In addition, the pressure and temperature corrections for 
binary diffusion coefficients are used [9]: 


Po T 1.5 
b= 2 (2) DÌ, (13) 


The molar fraction of nitrogen is calculated knowing that the 
molar fractions sum up to unity: 


Xn, = 1 — Xpo — Xo, (14) 


2.1.2. Cathode electrode 

The governing equations for the cathode electrode are the 
same as for cathode GDL with the exception that the conserva- 
tion equations of mass, species and charge have source terms. 
The mass and species equations have source terms because of the 
oxygen consumed and water produced in the fuel cell reactions, 
and charge equation because the protonic current is changed 
into electronic current. Thus, these conservation equations are 
rewritten as in Eqs. (15)-(19): 


Ce seo: P jMino (15) 

v- No, =- (16) 

ee (17) 
2F 

V + (>06 V pe,c) = je (18) 

V -(-omV¢m) = = je (19) 


The current production due to electrochemical reactions at 
the cathode is calculated from the Butler—Volmer equation: 


; . CO —a,F 
je = joc = exp ( n) (20) 
CO, 


RT 


The reference oxygen concentration is taken to be the case 
where only pure oxygen is present at the electrode, and thus 
the concentration term in the Butler—Volmer equation can be 
approximated as the molar fraction of oxygen, i.e.: 

CO, 


ref 
C Oz 


so (21) 


The cathode overpotential in Eq. (19) is defined as 
N = bec — m — Eo (22) 


2.1.3. Cathode boundary conditions 

The electronic potential of the cathode decreases at the inter- 
face between gas diffusion layer and electrode due to contact 
resistance. The potentials of electrode and gas diffusion layer 
are related to each other through the current density passing the 
interface by Ohm’s law giving a condition for both electrode and 
gas diffusion layer potentials at boundary X: 


OParLc _ g Pec _ Pec — PODL. 
dy e ay 


i= OGDL, y (23) 


Fcont,e 


The potential loss in the current collector is assumed to be 
negligible, and thus the only loss between the gas diffusion layer 
and current collector is due to contact resistance. Similarly as 
in Eq. (23) the contact resistance between the current collecting 
rib and GDL at boundary I yields: 


IPGDL,c _ PGDL,c — $0,c 


(24) 


i = —OGDL,) 
ð Fcont, gr 

No electronic current passes through the interfaces between 
the GDL and the channel (boundary II), and electrode and mem- 
brane (boundary XI), and thus 


IPGDL,c _ e.c 
dy ay 
The protonic potential is set continuous over boundary XI 

because it is assumed that there is no contact resistance for pro- 

tonic current between membrane and electrode. Gas diffusion 
layer does not conduct protons, and thus no protonic current 
passes through boundary X yielding: 


dpm 
dy 

It is assumed that the gas mixture is at ambient pressure in the 
channel and that the gas mixture is always ideally mixed, leading 
into fixed boundary conditions for pressure and species molar 
fractions at boundary II. Typically in 2D models, a standard value 
of 0.21 is used for molar fraction of oxygen. However, such a 
condition exists only at the very beginning of the inlet channel 
where no oxygen has been consumed, or when an infinite air 
stoichiometry is used. For this reason, it is assumed here that 
the molar fraction of oxygen corresponds to the average value 
in the middle of the channel when dry air with a stoichiometry 
of 2 is fed into the cell. This leads into fixed values of approx- 
imately 0.153 and 0.077 for oxygen and water molar fractions, 
respectively. The fixed molar fraction of nitrogen at boundary II 
is calculated from Eq. (14). 

There is no mass transfer at the interfaces between the gas 
diffusion layer and current collecting rib (boundary I), and elec- 
trode and membrane (boundary XI), and thus 


=0 (25) 


=0 (26) 


Ù- èy = (27) 
Ni-éy =0 (28) 


Because the anode side mass transfer was neglected, it is 
assumed for simplicity that also the water does not penetrate the 
membrane. 

Finally, symmetry boundary conditions are applied on the 
boundaries III and IV, i.e.: 


ù- e =0 (29) 
N;-é, =0 (30) 
O0@GDL Ë dpe c Ipm 
— = — = = 1 
əx Ox Ox n BL) 


2.1.4. Membrane and anode 
It was assumed that there is no mass transfer of water in the 
membrane. Thus, the only governing equation at the membrane 
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is the conservation of charge: 
Vom = 0 (32) 


Due to significantly faster electrode kinetics and smaller mass 
transfer limitations compared to the cathode, the conservation 
equations of mass, momentum, and species are not solved at 
the anode. Thus also at the anode gas diffusion layer the only 
equation to be solved is the conservation of charge: 


a a> 
fona, ) =, 
dy 


IPGDL,a = 
ex 
Ox 


V. ( OGDL,x OGDL,y (33) 


At the anode electrode the electronic current is consumed and 
protonic current produced yielding: 


V+ (-ceVoe,a) = —Ja (34) 
V - (—0mV m) = ja (35) 
where the current production at the anode is calculated from the 
Tafel equation: 

— Jo,az F 
a ie 


(be,a — m) (36) 


2.1.5. Membrane and anode boundary conditions 

It is assumed that there is no contact resistance for protonic 
current at the interfaces between membrane and electrode, and 
thus the protonic potential is continuous over boundary XII. 
Boundary conditions for electronic current at the anode elec- 
trode and GDL are similar to those at the cathode. The contact 
resistances at the interfaces between electrode and GDL, and 
current collecting rib and GDL yield conditions for boundaries 
XIII and VII: 


roy IPGDL,a T Ibe,a _ PGDL,a — e,a (37) 
Di dy j dy cont,e 
: O@GDL,a 90a — PGDL,a 
_ . a _ ®o, f 38 
E ð Fcont, gr (38) 


Because the mass transfer at the anode was neglected, the 
interface between the GDL and gas channel (boundary IX) has 
only a no current condition: 


OPGDL,a 
dy 


Finally, the symmetry boundary conditions apply again at 
boundaries V-VII: 


dpm Pea IPGDL,a 
əx ax æ 


=0 (39) 


=0 (40) 


2.1.6. Effect of inhomogeneous compression 

Due to the inhomogeneous compression of the gas diffusion 
layer its properties are changed. The affected GDL properties 
are porosity, permeability, in- and through-plane bulk conduc- 
tivities, and contact resistance between gas diffusion layer and 
electrode interface. These changes are taken into account in 
the modeled cathode domain B of Fig. 1b. Cathode domain 
A remains unchanged compared to the base case of Fig. la. 
Because mass transfer is considered negligible at the anode, the 


effect of the porosity and permeability changes can be neglected 
there. The only differences at the anode compared to the base 
case are the varying bulk conductivities of GDL and contact 
resistance at GDL/electrode interface due to inhomogeneous 
compression. 

The results of the measurements in the experimental part [6] 
imply that the gas diffusion layer is very little compressed in the 
middle of the channel, and that the total change from the original 
uncompressed volume remains small. Thus, the thickness of the 
gas diffusion layer in the cathode domain B is modeled with a 
logarithmic curve having a maximum that corresponds to 10 ym 
compression from the original thickness ho. The curve was fitted 
so that it coincides with the constant compressed thickness Acomp 
at the point where the gas channel ends and current collecting 
rib begins, and equals the maximum thickness in the middle of 
the channel, i.e. at the right boundary III. The resulting function 
for gas diffusion layer thickness implemented into the model 
geometry coordinates is 


h(x) [m] 


hcomp, xeA (41) 
~ | 19.30314 log((x — 0.0005) x 10°+ 1) x 1076 + hcomp, ¥€B 


It is assumed that the change in thickness under compression 
is due to change in volume of pores, not in volume of bulk 
material. Thus, the porosity of the GDL can be calculated from 
the thickness as 

h(x) — hmi 
(x) = eq 1) = min (42) 
ho — hmin 
where hmin equals the minimum thickness when there is only 
bulk material left, i.e.: 


hmin = (1 — €0)ho (43) 


A third degree polynomial fit was made with the least square 
sum method to the permeability data from the measurements 
[6], and the yielding function (fitting accuracy of R? = 0.994) 
for GDL is 


k(x) [m7] = —1.700 x 107!! + 2.760 x 1077h) 
— 1.484 x 107FA(x)? + 2.754h(x)? (44) 
The GDL in- and through-plane bulk conductivities were 


modeled as linear fits from the experimental data (fitting accu- 
racies of R? = 1.000 and 0.975, respectively), and they were 


OGDL, x(x) [Q7! m7!] = 6896 — 1.159 x 107A(x) (45) 
OGDL,y(x) [27 m~"] = 3285 — 8.385 x 10°A(x) (46) 


An exponential fit for the experimental data (R? = 0.983) of 
the contact resistance between GDL and current collecting rib 
gave 


Feont,gr(x) [2 m*] = 5.83 x 107!° exp(2.06 x 10*h(x)) (47) 


Because the experimental data of the contact resistance at 
GDL/electrode interface was unreliable, it was determined from 
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Table 1 

Dimensions of the modeled geometries 

Parameter Symbol Value 
Channel and rib width d 1 mm 
Uncompressed GDL thickness ho 380 um 
Compressed GDL thickness lcomp 250 um 
Electrode thickness 10 um 
Membrane thickness 25 um 


the contact resistance of GDL/rib interface. A correction factor 
for the Nafion content of the electrode, typically approximately 
30 vol.%, was used because Nafion does not conduct electrons. 
Thus, the used function for contact resistance was 

Feont,e(X) [2 m?] = 


Teont,gr(X) = 1.429rcont,gr(X) (48) 


1 
1— 0.3 

The above-listed parameter values for domain A and for the 
base case were calculated from the fitted equations in order to 
have continuous and similar parameter values between different 
domains and models. 


2.2. Parameters and model solving 


The dimensions of the modeled geometries are given in 
Table 1. The constants and parameters used in the model are 
listed in Table 2. Standard textbook values for constants and 
typical values found in the PEMFC modeling articles for fuel 
cell parameters are used when a reference is not given. 

The modeling was done using a commercial finite element 
method program COMSOL Multiphysics version 3.2b (formerly 
known as FEMLAB) with a parametric nonlinear direct (UMF- 
PACK) solver. When solving the model, the cell voltage was used 
as a fixed parameter by setting the potential of anode current col- 
lector to zero and the potential of cathode current collector to cell 


Table 2 

Constants and parameter values 

Parameter Symbol Value 

Ambient pressure Po 101,325 Pa 

Binary diffusion coefficient O2,H2O DÈ, mo 3.98 x 1075 m? s7! 
Binary diffusion coefficient O2,N2 Da i 2.95 x 1075 m? s7! 
Binary diffusion coefficient H20,N2 DÌ ok. 4.16 x 1075 m? s7! 
Conductivity of electrode Oe m 300 27! m7! [6] 
Exchange current density, cathode Joc 20 x 10° A m~? 
Exchange current density, anode Joa 1.7 x 10°? A m~? 
Faraday constant F 96,487 A s mol”! 
Gas constant R 8.314 J mol~! K7! 
Molar mass of oxygen Mo, 0.032 kg mol~! 
Molar mass of water My,0 0.018 kg mol! 
Molar mass of nitrogen My, 0.028 kg mol~! 
Permeability of electrode ke 1.26 x 107!3 m? [10] 
Porosity of uncompressed GDL £0 0.84 [11] 

Porosity of electrode Ee 0.4 

Protonic conductivity Om 5927! m! 
Reaction symmetry factor Oy 0.5 

Reversible cell potential Eo 1.23 V 
Temperature T 323.15K 

Viscosity of air H 1.9 x 107 kgm! s7! 


voltage. The used mesh consisted of 24,089 elements for base 
case and 24,420 elements for the case where inhomogeneous 
compression was taken into account. The respective degrees of 
freedom were 36,383 and 46,443. 


3. Results 


The polarization curves of the both simulated geometries are 
illustrated in Fig. 2. These were achieved by changing the cell 
voltage in steps of 0.1 V and calculating the average current 
density at each voltage over boundary X as 


; 1 T 1 
lave = Q 
a d Jo Fcont,e = 


There are no significant differences between the modeled 
cases at practical cell voltages implying that the overall cell 


PGDL,c) dx (49) 


1* T T T T T 
—— Base case 
0.9 —<— Inhom. comp. |] 


Voltage / V 


0.2 1 1 1 L 4 
0 300 600 900 1200 1500 
Current density / mA cm? 
Fig. 2. Polarization curves. 
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Fig. 3. Oxygen molar fraction at 0.4 V for the base case. 
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Fig. 4. Oxygen molar fraction at 0.4 V when the inhomogeneous compression 
is taken into account. 
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Fig. 5. Current density distribution and oxygen molar fraction at the 
GDL/electrode interface at 0.4 V for the base case. 
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Fig. 6. Current density distribution and oxygen molar fraction at the 
GDL/electrode interface at 0.4 V when the inhomogeneous compression is taken 
into account. 
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Fig. 7. Current production profiles at the electrode at 0.4 V for the base case. 
The profiles are current production rates in y-direction drawn at every 0.1 mm 
in x-direction. The arrow in the figure shows the direction of increasing x-axis. 
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Fig. 8. Current production profiles at the electrode at 0.4 V when the inhomo- 
geneous compression is taken into account. The profiles are current production 
rates in y-direction drawn at every 0.1 mm in x-direction. The arrows and labels 
in the figure show the direction of increasing x-axis and corresponding values 
in millimeters. 
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Fig. 9. Current density profile in the cathode electrode at 0.4 V when the inho- 
mogeneous compression is taken into account. 
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Fig. 10. Current density distributions at the GDL/electrode interface at 0.4V with varying contact resistance values between GDL and electrode when the 
inhomogeneous compression is taken into account: (a) conte = 0.01 rcontgr> 0.1 rcontgrs and 0.5/cont,gr3 (b) reont,e = Ycont,gr> 1/1 — 0.3)rcontgr> and 3Fcont,gr- 


performance is not significantly affected by the inhomogeneous 
compression of GDL. 

The differences in oxygen molar fractions between the mod- 
eled cases were also quite small. As an example, the contour 
plots of oxygen molar fraction at cell voltage of 0.4 V are illus- 
trated in Figs. 3 and 4. When the inhomogeneous compression is 
taken into account, the equimolar lines are slightly more vertical 
and shifted towards the parts below the current collecting rib. For 
example, there is only a 5% difference in oxygen molar fraction 
at the GDL/electrode interface below the middle of the channel 
at 0.4 V. There is also no significant difference when the liq- 
uid water saturation begins (at the electrode/membrane interface 
below the middle of the rib) between the modeled cases, being 
519 mV for the base case and 513 mV when the inhomogeneous 
compression is taken into account. 

The small differences in mass transfer are due to relatively 
high open porosity and permeability of the used GDL. With 
another GDL material that has a microporous layer and denser 
structure, the differences in mass transfer are most probably 
more pronounced. In addition, the capillary movement of liq- 
uid water was neglected, and thus the mass transfer at low cell 
voltages is somewhat distorted. 

Even though there were no significant differences in the over- 
all cell performance and molar fractions between the different 
modeled cases, the current density distribution is significantly 
affected by the inhomogeneous compression. The current den- 
sity profiles at GDL/electrode surface at 0.4 V are illustrated 
in Figs. 5 and 6. Also the oxygen molar fractions at the same 
interface are illustrated in the figures. 

The shape of the current density distribution follows quite 
much the oxygen molar fraction profile in the base case, which is 
a very typical modeling result. When the inhomogeneous com- 
pression of GDL is taken into account, there is a significant 
increase in the current density below the position where the 
channel begins. This increase is not due to numerical inaccura- 
cies of the solution, because the shape of the peak was unaffected 
by the density of the mesh. 

The current production profiles, illustrated in Figs. 7 and 8, 
show that there are significant differences in the reaction rates 
between the modeled cases. The reaction rate is decreased below 
the channel when the inhomogeneous compression is taken into 
account because of higher resistive losses. However, the differ- 


ences in reaction rates are quite moderate compared to current 
density distribution in order to explain the observed peak in 
Fig. 6. 

The reason for the peak is illustrated in Fig. 9, where the 
current density profile in the cathode electrode is plotted. A sig- 
nificant portion of the current produced in the parts below the 
channel flows in-plane in the electrode and enters the GDL from 
the part where the contact resistance is decreased. This phe- 
nomenon obeys the second Kirschoff law: the amount of current 
going through a certain route is inversely proportional to the total 
resistance of that route. Even though the contact resistance value 
between GDL and electrode was estimated from the contact 
resistance between GDL and graphite, the effect exists even if the 
value was highly overestimated. This is due to the fact that when 
the total resistive losses under the channel are increased, a bigger 
portion of the produced current flows laterally in the electrode 
towards the smaller resistance. This phenomenon is illustrated 
in Fig. 10, where the current density distribution is calculated 
with several different contact resistance values ranging from 
0.01 rcont.gr tO 37cont,gr- The current density distribution is some- 
what smoothened when the contact resistance is decreased, but 
even at two orders of magnitude smaller value, the distribution is 
still highly peaked. The average current density values between 
different cases vary from 966 to 1113 mA cm~? between the 
highest and lowest contact resistance values, respectively. 


4. Summary and discussion 


This paper focused on modeling the effects that the inho- 
mogeneous compression of gas diffusion layer has on the 
performance of a PEMFC. Model took into account the multi- 
component mass transfer in the cathode components and charge 
transfer in all of the cell components. Model was isothermal 
and the capillary movement of liquid water was not taken into 
account. The experimental parameters evaluated in Ref. [6] were 
used in the model, and the results were compared with a con- 
ventional model that excludes the effects of inhomogeneous 
compression. 

There were no significant differences in the overall cell per- 
formance between the modeled cases. In addition, the mass 
transfer was not significantly affected by the inhomogeneous 
compression, which was due to the used highly porous and per- 


T. Hottinen et al. / Journal of Power Sources 171 (2007) 113-121 121 


meable GDL. This may not be the case when a denser GDL 
material with a microporous layer is used, because then the 
differences in mass transfer are pronounced. In addition, the 
capillary movement of liquid water was neglected and thus the 
mass transfer at low cell voltages was somewhat distorted. 

The effect of inhomogeneous compression on the reaction 
rate was evident. The reaction rate was decreased under the 
channel because of higher total losses caused by increased 
bulk and contact resistances. Besides affecting the reaction rate, 
the effect of inhomogeneous compression on current distribu- 
tion was tremendous. The current density distribution on the 
GDL/electrode interface was peaked at the parts below the edge 
of the channel. This was due to redistribution of the current 
profile in the electrode. A significant portion of the current 
flowed in in-plane direction in the electrode, and entered the 
GDL below the rib where there was significantly lower contact 
resistance. This phenomenon was investigated with several dif- 
ferent contact resistance values between GDL and electrode, and 
even with very small values the current density distribution was 
significantly peaked. 

In the model it was assumed that there is a sharp edge in the 
shape of the GDL at the rib/channel interface. Even though the 
measurement results in Ref. [6] implied that the GDB is virtually 
not compressed below the channel, the edge is not necessarily 
that sharp in reality. This is especially the case when molded 
composite flow-field plates with slightly rounded corners are 
used. However, this causes that the changes in GDL properties 
under the channel are not that drastic only under the rounding and 
thus has only a small effect on the current density distribution 
by slightly widening the peak to the right and rounding the tip 
of the peak. 

The observed current density peak can have tremendous 
effects not only on current density distribution, but also to tem- 
perature distribution inside the cell. According to the analogy 
between charge and heat transfer, it can be assumed that also a 
significant portion of the heat produced in the electrode below 
the channel flows in in-plane direction. This means that there has 
to be a lateral temperature gradient within the electrode causing 
a possible hot spot below the channel. In addition, the Ohmic 
heating at the place where most of the current enters the GDL 
causes another possible hot spot below the place where the rib 


begins. This uneven temperature distribution can have signifi- 
cant effects on cell lifetime and the local cell performance. In 
the optimization of the PEMFC design the effect of these possi- 
ble hot spots should be minimized by minimizing the parts that 
have low compression pressure, but in a way that efficient mass 
transfer is simultaneously ensured. 

It was shown in this paper that the inhomogeneous compres- 
sion of GDL cannot be neglected. In addition, in order to reveal 
the possible hot spots caused by the lateral current and heat flow 
in the electrode, also the electrodes have to modeled as a separate 
domains, i.e. modeling the electrodes as boundary conditions as 
sometimes is done is not a valid approach. In order to achieve 
reliable estimates for the temperature distribution, the energy 
equations have to be included in the model. Before this can be 
accomplished, accurate evaluation of the GDL bulk and contact 
heat transfer parameters as a function of thickness is required. 
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